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ABSTRACT 


Renewed  interest  in  slip  flow  within  the  realm  of  continuum 
physics  has  resulted  in  a  numerical  study  on  generation 
and  propagation  of  vorticity  in  a  flow  past  thin  elliptic 
cylinders  under  the  perfect-slip  condition.  Computer  results 
show  that  surface  vorticity  is  more  concentrated  at  the  edges 
for  perfect  slip  than  for  nonslip,  due  to  the  surface 
curvature.  Flow  separation,  instability,  and  vortex  shedding 
occur  also  under  the  perfect-slip  condition.  Although  drag 
is  drastically  reduced  for  slender  body  configurations,  drag 
coefficients  for  a  plate  normal  to  the  flow  are  almost  the 
same  for  nonslip  and  perfect  slip,  at  least  for  Reynolds 
number  10.  All  computations  were  carried  out  for  Reynolds 
numbers  below  200. 
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REMARKS  ON  THE  GENERATION  OF  VORTICITY 


Current  experiments^  on  special  solid-surface  treatment  which 
may  prevent  fluids  from  adhering  to  the  surface  of  solid  bodies  have 
prompted  a  study  of  the  fluid  dynamic  implications  of  flows  under  the 
perfect-slip  condition.  This  type  of  motion  must  be  distinguished 
from  the  slip  flow  known  in  rarefied  gas  dynamics,  where  the  mean 
free  path  lengths  of  the  gas  molecules  are  considered  to  be  of  the 
order  of  the  body  dimension.  Here,  perfect-slip  flow  is  investigated 
within  the  framework  of  continuum  physics,  that  is,  the  mean  free 
path  length  is  very  small  with  respect  to  the  body  length.  Still, 
the  adherence  properties  must  be  explained  on  the  molecular  level, 
and  they  enter  the  flow  theory  of  a  continuum  only  as  boundary  con¬ 
ditions.  This  paper  does  not  discuss  the  cause  of  nonslip  and  does 
not  describe  methods  of  producing  slip  flow.  Rather  it  deals  with 
the  change  in  the  flow  behavior  when  the  nonslip  condition  is  replaced 
by  the  perfect-slip  condition.  The  solution  to  this  problem  is  not 
only  of  engineering  interest  but  furnishes  additional  insight  into 
the  generation  and  propagation  of  vorticity.  The  latter  aspect  is 
stressed  in  this  study. 

In  an  incompressible  homogeneous  fluid,  which  is  not  subjected 

to  nonconservative  forces,  vorticity  can  be  created  only  at  the 

boundaries  of  the  fluid.  For  a  body  moving  in  a  fluid  at  rest  at 

infinity,  vorticity  is  thus  produced  at  the  body  surface.  For 

simplicity  we  restrict  ourselves  to  plane  and  axisymmetric  flows 

and  use  the  intrinsic  coordinates  s  and  n,  where  n  denotes  the  arc  ' 

2 

length  along  streamlines, and  n  is  orthogonal  to  s,  n  =  n^  represents 
the  surface  of  the  body.  With  u  and  x  designating  the  velocity 
component  and  the  shear  stress  along  the  streamlines,  the  nonslip 
and  perfect-slip  conditions  at  the  surface  are  defined  by 


1.  References  are  listed  on  page  15 
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u  =  0, 

Nonslip 

(1) 

x  e  0,  y  £  0, 

Perfect-Slip 

(2) 

where  y  is  the  dynamic  viscosity. 

The  amount  of  vorticity  produced  at  the  surface  is  determined  by 
the  velocity  field.  Vorticity  defined  as  the  curl  of  the  velocity 
vector  is 

=  "  +  KU>  (3) 
The  shear  stress  in  terms  of  u  is 

t  -  -y  (|^  +  ku)  ,  (4) 

where  k  designates  the  curvature  of  the  surface.  In  plane  motions  o)! 
is  the  vorticity  component  normal  to  the  plane;  in  axisymmetric  flows 
03 1  is  the  azimuthal  component.  Under  the  nonslip  condition,  Equation 
(1),  vorticity  and  shear  stress  at  the  surface  are  u)ts=  x/p  =  -3u/3n, 
whereas  the  perfect-slip  condition,  Equation  (2),  requires  o)f  =  2ku  = 
-23u/3n.  Hence,  under  the  perfect-slip  condition  vorticity  can  only 
occur  when:  the  curvature-  k  is  nonzero.  ;; 

SLOW  MOTION  AROUND  A  SPHERE 

Vorticity  generation  and  diffusion  under  either  nonslip  or  perfect 

slip  may  be  illustrated  by  the  closed-form  solution  for  the  flow  past 

a  sphere  at  vanishing  Reynolds  number.  This  well-known  solution 

consists  of  two  parts;  one,  called  the  Stokeslet,  deals  with  the 

vorticity  distribution  and  the  other,  a  potential  flow,  satisfies  the 

boundary  conditions.  The  potential  flow  consists  of  a  dipole  and  a 

parallel  flow. 

.  U  .  2.  /A  ,  B  ,  2. 

ip  =  -  2  Sln  <f>  (Ar  +  —  +  r  ),  (5) 

oj  1  =  U  sin<j)  — r  ^  0*  (6) 

r 

Here,  ijj  is  the  stream  function  in  the  spherical  polar  coordinate 
system  (r,  <f>,  X),  and  U  is  the  uniform  flow  velocity  far  from  the  body. 
The  radius  of  the  sphere  is  unity.  The  constants  A  and  B  are  at  our 
disposal  to  fulfill  the  boundary  conditions  on  thesphere’s  surface. 
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A  remarkable  connection  exists  between  dipole  and  Stokeslet.  In 
a  dipole  vorticity  is  concentrated  at  the  singularity  r  =  0.  If  this 
source-sink  ceases  to  produce  vorticity  at  the  time  t'  =0,  the  dipole 
decays  according  to 


*D  *  '  UB 


t«f  (§ .-I,"*,, 


—  UBr  sin  <J>  (  vt')”5^2  e” ^  ^ , 
4/rT 


where  5  =  r/A>  t'  and  v  is  the  kinematic  viscosity.  Equation  (9) 
was  found  by  Phillips^,  Equation  (8)  by  the  author^.  If  the  dissipating 
dipole  is  integrated  over  time,  the  accumulation  of  vorticity  results 
in, 


(t  *  —  t)  dt 


U  *2  ,  r  -.  /!  2v  -  2 

-  Br  sm  <p  [  1-(1-  — )  erf  (-) - 

l  rlLr  r— 


[1-erf  (f)  +1  e-«2/4]. 

2  7T 


e-?2/4  ]  , 


These  expressions  yield  the  Stokeslet  for  tT  =  °°,  r  ^  0. 

The  nonslip  or  perfect-slip  condition  enters  when  the  constants 
A  and  B  in  Equations  (5)  and  (6)  are  determined.  These  constants  do 
not  change  the  characteristics  of  the  dipole  and  the  Stokeslet.  We 


obtain  for 
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Nonslip:  A  =  -  — ,  B  =  —  (Stokes,  1851), 

Perfect  Slip:  A  =  -1,  B  =  0  (Basset,  1888). 


The  drag  for  the  sphere  of  unit  radius  is  6ttijU  under  the  nonslip 
condition,  4TryU  under  the  perfect-slip  condition^. 

PROPAGATION  OF  VORTICITY  UNDER  PERFECT  SLIP 

To  study  the  propagation  of  vorticity  by  means  of  diffusion  as  well 
as  convection  under  the  perfect-slip  condition,  solutions  of  the 


Navier-Stokes  equations  are  constructed  numerically/  for  nonzero  Reynolds 
numbers.  The  computer  program  developed  to  perform  the  calculations 
is  described  in  Reference  6.  Sufficient  notations  and  remarks  are  given 
in  the  following  text  such  that  the  reader  can  understand  the  results 
which  are  presented  in  Figures  2  through  8. 

The  computer  program  is  designed  to  calculate  two-dimensional, 
parallel,  unsteady  viscous  flows  of  incompressible  fluids  past  elliptic 
cylinders  under  various  angles  of  attack.  Potential  flow  is  selected 
as  the  initial  condition,  thus  simulating  the  abrupt  start  of  the  body. 
With  time  either  a  steady  state  or  a  periodic  motion  develops  (for 
the  range  of  Reynolds  numbers  under  consideration).  The  infinite  flow 
field  is  represented  by  a  finite  grid  system  which  is  based  on  the 
elliptic  coordinates  n  and  0  defined  by 

x  +  iy  =  a  cosh  (n  +  i0),  a  >  0,  i  =  -1.  (13) 

The  focal  distance  is  denoted  by  a,  and  the  body  contour  is  described 
by  n  =  =  0.1.  Figure  1.  The  space  increments  are  An  =  0.05, 

A0=  ir/30  or  tt/40,  and  the  number  of  mesh  points  is  75  x  60  =  4500  or 

75  x  80  =  6000.  The  Reynolds  number  is  defined  by  Re  =  dU/v,  where 
d  =  2a  cosh  n^.  The  time  t'  is  made  dimensionless  by  t'  =  ta/U.  The 
computer  programs  for  nonslip  and  perfect-slip  flows  differ  in  the 
calculation  of  the  surface  vorticity,  surface  velocity,  pressure,  and 
the  drag,  lift,  and  moment  coefficients. 

The  following  cases  have  been  investigated  and  compared  with  the 
corresponding  solutions  for  nonslip  flows:  (a  =  angle  of  attack) 

a  =  0°,  Re  =  200,  symmetric  flow,  almost  steady  state,  75  x  60  grid 

a  =  90°,  Re  =  .10,  symmetric  flow,  almost  steady  state,  75  x  60  grid 

a  =  45°,  Re  =  200,  asymmetric  flow,  development  of  vortex  street, 

75  x  80  grid 
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Figure  1  -  Elliptic  Coordinate  System  and  Definition  of  Angle  of  Attack 

Surface  vorticity  and  surface  pressure  along  the  body  contour  for 
symmetric  flows  are  displayed  in  Figures  2  through  5.  Vorticity 
is  more  concentrated  at  the  edges  for  perfect  slip  than  for  nonslip* 
Away  from  the  edges  vorticity  is  almost  zero  for  perfect  slip.  This 
reflects  the  fact  that  the  occurence  of  vorticity  depends  on 
the  surface  curvature.  In  the  case  of  the  thin  elliptic  cylinder 
situated  perpendicular  to  the  flow,  the  surface  vorticity  changes  its 
sign  immediately  behind  the  edges,  thus  indicating  flow  separation. 

The  streamline  patterns  which  are  presented  in  Reference  7  verify  this 
and  exhibit  wakes  with  twin  vortices  which  are  remarkably  similar  for 
nonslip  and  perfect  slip.  The  fact  that  flow  separation  occurs  for 
perfect-slip  flow  is  surprising.  Obviously,  the  vanishing  of  the 
shear  stress  at  the  surface  does  not  prevent  flow  separation  and  can¬ 
not  be  used  as  a  criterion  for  flow  separation  under  the  perfect-slip 
condition.  However,  the  definition  for  flow  separation, 

U>1  =  °,  <ff  )1<  0  (14) 

8 

based  on  the  vorticity  holds  for  both  nonslip  and  perfect  slip. 

The  peaks  in  the  curves  of  the  surface  pressure  p^  at  the  edges 

are  also  more  pronounced  for  perfect  slip  than  for  nonslip.  Both  to- 

2  1 

and  p^  are  used  to  compute  the  coefficients  =  Drag/  (p/2)U  d, 

CL  =  -  Lift/ (p/2)U2d,  and  CM  -  -  Torque/ (p/2)U2d2.  With  CD  =  CDp  + 
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Figure  2  -  Surface  Vorticity  under  Nonslip  and  Perfect  Slip  for 
Re  =  200,  a  =  0°,  n,  -  0.1  at  almost  Steady  State 
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Figure  3  -  Surface  Pressure  under  Nonslip  and  Perfect  Slip 
for  Re  =  200,  a  =  0°,  =  0.1  at  almost  Steady  State 
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Figure  4  -  Surface  Vorticity  under  Nonslip  and  Perfect  Slip 
for  Re  =  10,  ot  =  90  ,  “  0.1  at  almost  Steady  State 
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where  and  are  the  drag  coefficients  due  to  pressure  and  friction, 
Dr  DJc 

respectively,  the  results  are: 


a  =  0°,  Re  =  200,  ri ^  =  0.1:  Nonslip 

Perfect-Slip: 


a  =90°,  Re  =  10  ,  =  0.1,  t  =  12 


Nonslip : 
Perfect  Slip: 


CD  =  0.55,  CDP  =  0.07 
CD  =  CDP  =  °’°°92 


CD  =  6.00,  CDp  =  5.55 
CD  -  CDp  =  5.73 


Hence,  the  drag  coefficient  for  the  thin  elliptic  cylinder  parallel/’ to 
the  flow  decreases  60  times  when  the  nonslip  condition  is  replaced 
by  the  perfect-slip  condition.  However,  and  this  is  really  surprising, 
for  the  elliptic  cylinder  normal  to  the  flow,  the  coefficients  are 
almost  the  same.  The  reason  is  that  perfect  slip  reduces  drag  con¬ 
siderably  only  for  bodies  with  dominating  wall-shear  stress.  Flow 
displacement,  on  the  other  hand,  is  scarcely  affected  by  the  change 
in  the  surface  condition. 

Once  vorticity  is  produced,  it  propagates  by  diffusion  and  con¬ 
vection,  independently  of  how  it  is  produced.  Only  in  the  immediate 
neighborhood  of  the  surface  is  propagation  directly  influenced  by 
the  surface  condition.  Taylor-series  expansions  around  a  surface 
point  reveal  that  only  diffusion  contributes  to  the  leading  terms 
under  the  nonslip  condition,  whereas  both  diffusion  and  convection 
participate  when  perfect  slip  is  assumed. 

The  case  a  =  45°,  Re  =  200  illustrates  the  behavior  of  unsteady 

asymmetric  flows  under  perfect  slip.  Figures  6  and  7  show  patterns 
of  streamlines  and  lines  of  constant  vorticity  for  the  initial  phase 
of  a  developing  Karman  vortex  street.  The  pictures  are  remarkably 
similar  for  both  nonslip  and  perfect  slip  except  in  the  central  region 
of  the  body  surface.  Under  perfect  slip  an  initial  vortex  is  generated 
as  in  the  nonslip  flow  giving  rise  to  a  flow  circulatipn  (in  the  sense 


11 


NONSLIP 


PERFECT  SLIP 


Figure  7  -  Sequence  of  Lines  of  Constant  Vorticity  Corresponding  to 
the  Streamline  Patterns  in  Figure  6 


of  the  Lanchester-Prandtl  hypothesis  for  potential  flow).  In  Figure  8 
Cp,  C^>  and  are  plotted  versus  t.  The  orders  of  magnitude  for 
nonslip  and  perfect  slip  are  the  same.  Also  the  Strouhal  numbers^ 
agree  well.  A  difference  occurs  immediately  after  the  start.  Whereas 
for  nonslip  the  sudden  start  causes  infinite  values  of  and  and 
a  zero  value  of  C^,  the  situation  is  reversed  for  perfect  slip. 


Figure  8  -  Drag,  Lift,  and  Moment  Coefficients  Versus  Time  for  Nonslip 
and  Perfect-Slip  Flow  with  Re  =  200,  a  =  45°,  n^=  0.1 
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